A new ripplon branch in He II 
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We analyse the dispersion relation of ripplons, on the surface of superfiuid helium, using the 
dispersive hydrodynamics approach and find a new ripplon branch. We obtain analytical equation 
for the dispersion relation and analytic expressions for the limiting cases. We discuss where ripplons 
can exist in the energy-wavenumber plane. A numerical solution for the ripplon dispersion curve 
is obtained in the allowed regions. The new ripplon branch is found at energies just below the 
instability point. 



I. INTRODUCTION 

Ripplons are quantised capillary waves on the free sur- 
face of superfiuid 4 He. At low frequencies the dispersion 
law gives a good way of measuring the surface tension of 
liquid helium [1, 2]. The temperature dependence of the 
surface tension is due to ripplons [ ] . Ripplons can be de- 
tected by neutron scattering, in a similar way to the bulk 
modes, and have been shown to exist up to wavenumber 
of 1.5 A , where the ripplon energy is equal to that of 
the roton minimum, 8.6 K, [4, 5]. Ripplons are the dom- 
inant scatterer of surface state electrons on liquid helium 
(see [G, 7] and references therein). Ripplons play a signif- 
icant role in the condensation [8] , evaporation and reflec- 
tion of atoms from liquid 4 He, [9, 10]. It has been sug- 
gested that ripplons are the most favourable excitations 
for the simulation of general-relativistic effects related to 
horizons of white holes [11]. 

Ripplons are well defined excitations and so have a 
reasonably long lifetime, but they can be scattered by 
phonons and rotons, and by other ripplons. At low ener- 
gies there can be three ripplon scattering, in an analagous 
process to three phonon scattering. However for low en- 
ergy ripplons at low temperatures, the ripplon lifetime 
is dominated to phonon scattering rather than ripplon- 
ripplon scattering [12]. 

In this paper we present a theoretical model of rip- 
plons which is good enough to account for the measured 
dispersion curve w(fc) but is simple enough to expose the 
underlying physics. It explains why ripplons only exist 
in certain parts of the (oj, k) plane, and why the disper- 
sion curve approaches the line uj = A rot at the top of 
a parabola, where A rot is the energy of the roton mini- 
mum. Moreover the model predicts a new roton branch 
at energies <^ 2A rot . 
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The theoretical model is a non-local hydrodynamical 
theory [13], developed by us in Ref. [14]. The actual 
physical characteristics of the liquid are introduced into 
the model by using the measured dispersion curve for 
bulk excitations. The justification of application of non- 
local hydrodynamics to the description of a quantum 
fluid at interatomic distances is given in [14], and in more 
detail in [15]. It is based on the fact, that in a quan- 
tum fluid the atoms are delocalized, as their thermal de 
Broglie wavelength is greater than the atom spacing. The 
general idea of description of a quantum fluid in terms 
of hydrodynamic variables at interatomic scales is being 
widely used. One of the first to exploit it was Atkins 
in 1959 [1G], when he introduced bubbles and snowballs 
of microscopic size in order to describe the mobility of 
electrons and ions in He II by means of methods of the- 
ory of continuous medium. The idea is utilized now in 
the variations of density functional theory (see for ex- 
ample [17]), and is implicitly assumed in other fields of 
research related to superfluidity. In particular, the vari- 
ables of continuous media are used lately for description 
of vortices in helium with dimensions of the order of the 
interatomic distances [18]. 

The quantum fluid, which is considered continuous at 
any length scales, can be described in terms of the vari- 
ables of continuous media, which satisfy the mass and 
momentum conservation laws. The equations of ideal 
liquid, which follow from them, do not form a closed 
system, and are supplemented by the equation of state, 
which specifies the functional relation between pressure 
and density. As shown in [14, 15], the equation of state 
for small deviations of the system from equilibrium in 
the general case is non-local, with some difference kernel 
h(r). 

The closed linear system of equations is brought to a 
non-local integro-differential wave equation with regard 
to pressure. The dispersion relation of the fluid f2(fc) then 
is determined by the Fourier transform of the kernel h(r) 
and can contain arbitrary degrees of k. In this paper we 
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start from the dispersion relation that approximates the 
experimental data for the spectrum of superfluid helium. 
The chosen f2(fc) gives explicit expression for the kernel 
h(r), which determines the equation of state and the non- 
local wave equation. 

In order to derive the dispersion relation of ripplons, 
we look for the solution of the nonlocal wave equation 
(1) in half-space, equipped with usual boundary condi- 
tions (2). The use of simple model of superfluid helium 
allows us to fully solve the problem of ripplons' spec- 
trum. The ripplons' dispersion equation is obtained in 
algebraic form, and its analytic solutions are derived in 
limiting cases. The analytic and numerical results of the 
paper in the region below the roton gap are consistent 
with the experiments, numerical computations and qual- 
itative estimates of other authors [4], [17]. 

A preliminary and partial account of this analysis is 
given in [ ] , where a very approximate function for the 
dispersion curve for bulk phonons and rotons was used. 
This was useful to show that the method had promise 
but it could not show the new ripplon branch or other 
features. Here we use a much more detailed function for 
the phonon-roton dispersion f2(fc), which gives an excel- 
lent approximation of the experimentally measured curve 
over the whole wave vector range. This allows us to de- 
scribe the ripplons in detail in the whole energy interval. 
This range includes the energy of the Pitaevskii instabil- 
ity at twice the roton minimum energy [20]. It is just 
below this energy where we find a new ripplon branch. 



II. EQUATIONS AND BOUNDARY 
CONDITIONS 



Let us consider the half-space z > filled by superfluid 
helium. Following the approach in Ref. [14], the liquid 
obeys the ordinary linearized equations of an ideal liquid, 
but the relation between the deviations P and p of pres- 
sure and density from the respective equilibrium values 
is nonlocal, with some difference kernel h(r) [14, 15]. 



AP(r,i)= / ^rifcflr-riOPtn.t), 26(0,00). (1) 



zi>0 

Here the dots denote derivatives by time. 

We assume that the interface is sharp enough so that 
the kernel h(r) is the same near the interface as it is in 
the bulk medium. The kernel is related to the dispersion 
relation of the bulk excitations ft(k) through its Fourier 
transform (see [14]) h{k) = k 2 /fl 2 (k). 

The dispersion relation fi(fc), is for the continuous 
medium which fills the half-space z > 0. It is nonlinear 
and has the distinctive form of the dispersion of phonons 
and rotons in superfluid helium (see Fig.l), and the char- 
acteristic length of the kernel is of the order of average 
interatomic separation. 

The equation (1) is supplemented by the boundary 
condition at the free surface. The pressure at the surface 



with surface tension a, must include the Laplace term, 
which for small deviations of the surface from equilibrium 
is 



P 



\dx 2 



dy 2 



(2) 



where £ is the z-coordinate of the points of the surface 
(see for example [13] or [21]). 

All the variables have the form ~ exp(ikr — iujt). For 
a given frequency tu and the projection of wave vector k 
on the plane of the free surface (x,y), k r , this expression 
becomes P = — erfc 2 £. The z-component of velocity of 
the surface is v z = ~iuj^. We can therefore write the 
boundary condition in the form 



ak 2 



P\ 



(3) 



III. EQUATION FOR THE DISPERSION 
RELATION OF RIPPLONS AND ANALYTICAL 
SOLUTIONS 

The dispersion relation, for undamped waves, can in 
general be written u> 2 = fl 2 (k 2 ). We express the function 
il 2 (k 2 ) as a polynomial of degree S in powers of k 2 , such 
that the only real zero of fl 2 (k 2 ) is at k 2 — 0. Near this 
point O 2 ~ k 2 . After some analysis, the Fourier image 
by r and t of the solution of equation (1) can be written 
in the form (see [15] and references cited there): 



P(k z ;k T ,uj) = 



C OM t(tJ,k T ) t-t, k z -k iz (0,1s. T ) 
kg- k lz (u>, kr) k ±± c k z -k iz (uj,k T )' 



(4) 



Here k = k T +e z k z is the wave vector. The product is 
taken over all the roots k z = ki Z of equation VL 2 {k 2 = 
k 2 +k 2 )—u 2 in the upper half-plane C + of the complex 
variable k z . The real roots are assumed to be shifted up 
from the real axis. The root ki z ((jJ,\t T ) is the phonon 
root, i.e. the one in C+, which continuously goes to zero 
at uj = and k T = 0. The prime superscript on the product 
designates that k\ z is excluded from the product. The 
function C ou t(a;,k T ) is the amplitude of the wave. 

Ripplons have solutions which are damped with dis- 
tance from the liquid's surface. These occur for all the 
kiz, that are not real. For a given value of w, the imag- 
inary ki z occur when k 2 is greater than any of the real 
roots k 2 (co) of equation f2 2 (fc 2 ) = u 2 . 

The dispersion relation for superfluid helium is shown 
in Fig.l. We assign subscripts to the positive roots of 
equation il 2 (k 2 ) =uj 2 in the ascending order of their ab- 
solute values, i.e. fci, &2 and k% are phonons, FIT and R + 
rotons respectively. 

The inverse Fourier transform of (4) is calculated from 
the residues in C+, and velocity is obtained from the re- 
lation dv/dt = —\7P/po where po is the equilibrium den- 
sity. Then the values of P and v z , on the surface z = 0, 
can be expressed in terms of the residues of P(k z ) and 
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k z P(k z ) at infinity. These are found directly by expand- 
ing P from (4). Then we obtain 



where 



— \k lz +^T{k lz ~k lz (Lo = 0)}\p\ z=0 . (5) 



i=1 



Comparing (5) with the boundary condition (3), we fi- 
nally obtain the equation for the dispersion relation u>(k T ) 
of ripplons: 

uj 2 = — ■ J ki z (u},k T ) + y^ [k iz (u),k T )-k iz (0,k T )] \ ■ 

p° 1 I U J 

This equation is written in terms of ki Z (uj,k T ) = 
\jk\ (iS) — k% £ C+, where kf(co) in turn are the roots 
of the polynomial equation Q 2 (k 2 ) = u 2 with real coeffi- 
cients. Then taking into account the condition that all 
the ki z are non-real, we can show, that the right part 
of the equation (6) is in fact real. Indeed, let us con- 
sider the equation fl 2 (k 2 = k 2 + k 2 ) — ui 2 with regard to 
k 2 . This equation is also a polynomial equation with real 
coefficients. So its roots with regard to k 2 are either neg- 
ative (no positive roots exist or there would be real ki z ) 
or break up into complex-conjugate pairs. The negative 
roots give imaginary ki Z and complex-conjugate pairs 
kf z ,k 2 z give pairs ki Z ,kj z in C+, such that ki z — —k* z , 
and for such a pair k iz + kj z = 2Hmk iz . Therefore the 
expression inside of the braces in (6) is imaginary and 
the right hand part is real. 

Further insight is obtained if we search in Eq. (6) for 
to 2 in the form of a series expansion in k T . On substi- 
tuting the expansions of all quantities by small k T and u) 
into the equation (6), we obtain within O(fef) 



Pa 



1- 



a k T 



2p s 2 8p 2 



Po 



ki 



Pqs z 



(7) 

where P= ^(0,0)^ • 

U) — 

In the limit of small k T and u, we replace the braces in 
Eq. (7) by unity, and so obtain the well-known disper- 
sion of capillary waves. The second and third terms in 
the braces, take into account compressibility. The fourth 
term, proportional to /3, expresses the influence of the 
(S — 1) roots ki with i > 1, which include the rotons. 
The latter give only small correction at small w, but 
with increasing frequency they become of the order of 
the phonon terms and, as it will be shown below, they 
determine the asymptotic behavior of the curve u>(k T ) in 
the proximity of uj = A rot . 

Of particular interest is the behaviour of the dispersion 
curve close to the roton gap u> = A rot . In this region the 
functions k2,3 Z (u,k T ) contain terms ~ \/ A rot — to, which 
leads to the same square-root singularity in the equation 
(6): (k 2z + fc3 Z )~ y/A rot — uj. Therefore it can be shown 
that the asymptote of oj(k r ) at uj~ A rot — has the form 



'po K 



2f,/(k 2 ot k 2 ); 



-ic = k lz (A rotl k c )+J2 k iz (A rot , i 

i=4 



•£*r,(o>*c); 

i=2 



(8) 



k rot is the momentum of a roton with energy A ro t, P is 
the "roton mass" and k c is determined by the condition 
uj(k T = k c ) = Arot- The numerical solution below, shows 
that this point indeed exists, is unique and k c <k rotl so 
a, b, c are all real. We see that the curve u)(k T ) approaches 
the value lo — A rot at the top of a parabola, and ends in 
the adhesion point, with zero derivative. There is no 
dispersion curve below A rot with k T > k c . In the theory 
of Ref. [ ] the asymptote was obtained qualitatively 
as one of the possible variants, from general quantum- 
mechanical considerations. 

In Fig.l we show the dispersion curve for bulk exci- 
tations. The areas bounded by the dispersion curve and 
the lines at constant energy at A rot , A max , and 2A rot are 
labelled A to G. We consider the helium surface to have 
a sinusoidal perturbation with wavenumber k T and then 
see if a stable ripplon solution is possible in each of the 
areas. A ripplon solution requires k\ z to be imaginary. 
In region A, k T < fei and k\ z = k\ — k 2 so k\ z is real, and 
for k\z real, Eq. (6) gives k T complex which means the 
ripplon decays, so no ripplon can occur in this area. The 
imposed perturbation would decay into phonons. In re- 
gion B, k\ z is not real and therefore Eq. (6) gives k T real, 
so ripplons can exist. They will lie on a ripplon disper- 
sion curve in this area. In region C, ki z k 2z and k 3z are 
real and the perturbation decays into phonons, R~ rotons 
and R + rotons. In region D, k 2z and k 3z are real and the 
perturbation decays into R~ rotons and R + rotons. In 
region E, k$ z is real and the perturbation decays into R + 
rotons. So there are no ripplons in these regions. In re- 
gion F there are no solutions to Eq. (6) because of the 
high gradient of the dispersion curve, so there are no rip- 
plons. In region G, the gradient is smaller and there are 
solutions to Eq. (6) so we predict that there are ripplons 
in this area. 



IV. NUMERICAL SOLUTION AND THE NEW 
RIPPLON BRANCH 

We find the ripplon dispersion curve numerically in the 
range of k T £ (0, 3) A" 1 . It is shown on Fig.l. We see that, 
as predicted, the computed curve at small wave vectors 

3/2 

is close to the classical k T law, but deviates from it at 
larger k T and approaches the level of A rot at the top of 
an inverted parabola at fc c = 1.27A -1 . 

There is no ripplon solution in the region to the right 
(see Fig.l) of the R + roton dispersion curve for u E 
(A rot , A rnax ). However, if we further increase uj while 
moving along the curve k^(uj), as we get close to the re- 
gion where Q(k)/k reaches its maximum, the curve turns 
down and aims to the point of instability [ ] in almost 
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FIG. 1. The dots are experimental data for the dispersion 
relation of bulk excitations Q(k) [23], and the thin line 
shows its analytic approximation that is used, with 5 = 18. 
The thick line shows the results of numerical solution of 
Eq. (6) for the ripplon dispersion curve oj(k T ). Large 
dots with error bars are experimental data for the ripplon 
dispersion [ ]. Two black dots in the high-energy part of 
the bulk spectrum designate the end points of the very 
high energy ripplon solution, and the insert graph shows 
its behaviour schematically, where 7 = 0.7- 10 _3 A _1 and 
6= 1.6 • 10~ 3 K . The different regions of the (u),k) plane 
are labelled A to G. The usual ripplons are in region B 
and the new ripplons are in region G. In the other regions, 
ripplons cannot exist, as described in the text 



a straight line, the left hand part of Eq. (6) starts to 
change slowly. At the same time, some of the functions 
kf(ui) for i > 3, which were slowly changing functions 
below A maxi begin to change fast and give significant 
contribution to right-hand part of Eq. (6), so here the 
ripplon solution is determined by the functions ki Z . 

Therefore, when we search for the solutions on the R + 
roton curve, we find two common points above A max . 
The first point is at 2.52 A" 1 and 16.65 K, which is almost 
exactly at the maximum of f2(fc)/fc, and the second point 
is 2.75 A -1 and 17.30 K. The ripplon dispersion curve 
w(£) T ) between these points lies near, but below, the bulk 
dispersion. The two common points are adhesion points. 
This is partly the reason for the deviation between the 
curves being extremely small. At k T = 2.6A~ 1 , midway 
between the two end points, the deviation is 1.6 • 10~ 3 K 
or 0.7 • 10~ 3 A -1 , which is too small to see on the scale 
of the main graph, so in the inset to Fig.l we show this 
region expanded. Nevertheless the separation of the rip- 
plon dispersion curve from the roton dispersion curve, is 
greatly exaggerated. 

We now show that the common points of the ripplon 
dispersion curve w(fc T ) and the R + roton dispersion curve 
k 3 (io) are adhesion points. The Eq. (6) can be rewritten 
as F(uj, k T ) = 0. There, the term k$ z in F(oj, k T ), goes to 
zero as yJk T — k^(uj) near the R + roton branch. So the 
gradient of F(ui,k T ) in the plane (uj,k T ) tends to infinity 
on the curve k T = k^^ui) and is directed normal to the 



curve. From the other side, the gradient of F(u>,k T ) is 
directed normal to the curve uj(k T ), which is its level 
curve F = 0. Therefore at the common points of the 
two curves, the angle between them is zero, and so these 
points are adhesion points. 

We now consider the penetration depth of the ripplon 
solution on this branch of the ripplon dispersion curve. 
At k T = 2.6 At 1 the value of k^ z is close to its maximum 
0.06 A -1 . Hence the penetration depth of the ripplon 
solution, determined by the term containing ^k% z , has 
a minimum value S ~ | fc 3 ^ | 1 ~ 16A, and at the end 
points it tends to infinity. The deviation between the 
two curves is second order in the small parameter k$ z 
because \k% z \ ~ \/k T — k^. 

The relatively large penetration depth of the solution 
means that macroscopic films or the surface of bulk he- 
lium are needed to observe these ripplons. They will not 
be seen on films of a few mononlayers. The high-energy 
ripplons should exist and be observable on saturated films 
of He II, which have a typical thicknesses of 300A. In the 
same way the penetration depth of the ripplons in region 
B (see Fig.l) tends to infinity when u tends to A rot . 
So the ending point of the dispersion curve (A rot ,k c ) 
can also be observed only in thick enough films of He II. 
Thus in the two most interesting regions of their disper- 
sion curve, the ripplon solutions have large penetration 
depth, much larger than the characteristic distances, of 
several monolayers thickness, of the changes of the den- 
sity profile at the free surface. The major part of the en- 
ergy of the wave is stored outside of the transition layer, 
it can be neglected, and the surface can be considered 
sharp (i.e. the kernel h(r) in Eq. (1) is the same as in 
the bulk fluid). 

The numerical calculations were carried out by repre- 
senting the function 2 (fc 2 ) as a polynomial of degree S, 
with 5=18 and 5 = 21. The high-energy ripplons solu- 
tion exists in both cases, and the end points do not differ 
substantially. The deviation between the curves k 3 (uj) 
and oj(k T ) remains very small, and is much less than the 
difference between the two polynomials. 

We have also checked for the stability of the solution 
with regard to changes in surface tension a, as at high 
curvatures of the surface its curvature dependence might 
be significant. If we decrease a, the end points of the 
ripplon dispersion curve slide apart along k^{uj). At the 
value a = 1.75 N/m (a = 3. 544 N/m at zero temperature 
[23]), the lower adhesion point reaches the point 2.44A -1 , 
16 K. If we increase a, then the end points slide towards 
each other and the solution disappears at a = 8.7 N/m. 
This shows that any realistic dependence a(k T ) can be 
inserted in Eq. (6) without complicating the numerical 
solution. 



V. CONCLUSION 

We have developed a non-local hydrodynamical model 
of ripplons on the free surface of superfluid 4 He, which 
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uses the measured dispersion curve for phonons and ro- 
tons to include the relevant characteristics of the liquid. 
The model gives a very good description of the measured 
ripplon dispersion curve, and it accounts for the parabolic 
approach of the dispersion curve to the line U) — A rot . 
The model gives the correct classical dispersion of capil- 
lary waves at long wavelengths. These good agreements 
validate the model. Furthermore the model is transpar- 
ent enough to expose the physical reasons why the low 
energy ripplon branch lies below the phonon dispersion 
curve and ends at the energy A rot . It is explained, with 
the aid of Fig.l, that surface perturbations, on the low 
k side of the phonon-roton dispersion curve, just decay 
to phonons or rotons and so stable ripplons can only be 
found on the high k side of the phonon-roton dispersion 



curve. For energies > A max , and on the high k side of 
the phonon-roton dispersion curve, ripplon solutions of 
Eq.(6) are possible at energies just below 2A rot . It can 
be shown that these ripplons are stable against decay into 
two ripplons with energy < A rot . That such ripplons so- 
lutions exist is unexpected and is a striking prediction 
of this analysis. We hope this prediction stimulates new 
experiments to detect these ripplons. 
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